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Global frequentist fits to the CMSSM and NUHMl using the MasterCode framework 
are updated to include the public results of searches for supersymmetric signals using 
~ 1/fb of LHC data recorded by ATLAS and CMS and - 0.3/fb of data recorded by 
LHCb. We also include the constraints imposed by the electroweak precision and B- 
physics observables, the cosmological dark matter density and the XENONIOO search 
for spin-independent dark matter scattering. Finally we also investigate the impact of 
a possible Higgs signal around 125 GeV. The various constraints set new bounds on 
the parameter space of the CMSSM and the NUHMl. We discuss the impact of the 
new results from SUSY and Higgs searches for these bounds and analyze the impact 
for a future liner e^e~ collider. 



1 Introduction 

The MasterCode collaboration pi-[TU] has reported the results of global fits to pre-LHC 
data [TH5] (see Ref. [TT] for other works) as well as including LHC 2010 data [5HTU] (see 
Ref. |12| for other works) in the frameworks of simplified variants of the minimal super- 
symmetric extension of the Standard Model (MSSM) [13] with universal supersymmetry- 
breaking mass parameters at the GUT scale. We consider a class of models in which R-parity 
is conserved and the lightest supersymmetric particle (LSP), assumed to be the lightest neu- 
tralino [13] , provides the cosmological cold dark matter [15] . The specific models studied 
have included the constrained MSSM (CMSSM) with parameters mo, TOi/2 and Aq denoting 
common scalar, fermionic and trilinear soft supersymmetry-breaking parameters at the GUT 
scale, and tan /? denoting the ratio of the two vacuum expectation values of the two Higgs 
fields. Other models studied include a model in which common supersymmetry-breaking 
contributions to the Higgs masses are allowed to be non-universal (the NUHMl), a very 
constrained model in which trilinear and bilinear soft supersymmetry-breaking parameters 
are related (the VCMSSM), and minimal supergravity (mSUGRA) in which the gravitino 
mass is required to be the same as the universal soft supersymmetry-breaking scalar mass 
before renormalization (sec Ref. [8] for an extensive list of references for those models) 

The impressive increase in the accumulated LHC luminosity combined with the rapid 
analyses of LHC data by the ATLAS |l6l[l7|, CMS [EHSO] and LHCb Collaborations [21] 
have been included in the analysis presented in Ref. [8]. Most recently, the ATLAS and 
CMS Collaborations have presented preliminary updates of their results for the search for 
a SM-like Higgs boson with ^ 5/fb of data [22]. These results may be compatible with a 
SM-like Higgs boson around Mh — 125 GeV, and they have been included in the analysis 
presented in Ref. [l]. 

Here we review the results of Refs. [8l[9], which are focused on the CMSSM and the 
NUHMl. We discuss their impact on the physcis prospects for a future linear e^e" collider, 
in particular the ILC, with a center-of-mass energy around 1 TeV (the ILC(IOOO)). 
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2 The MasterCode framework 

We define a global likelihood function, which combines all theoretical predictions with 
experimental constraints (except the latest LHC SUSY and Higgs searches): 



Here N is the number of observables studied, Ci represents an experimentally measured 
value (constraint) and each Pi defines a prediction for the corresponding constraint that 
depends on the supersymmetric parameters. The experimental uncertainty, tT(Ci), of each 
measurement is taken to be both statistically and systematically independent of the cor- 
responding theoretical uncertainty, cr(Pi), in its prediction. We denote by x^(-^?i) a-nd 
X^(BR(i?s — >■ the contributions from the two measurements for which only one- 

sided bounds were previously included. 

We stress that the three SM parameters /sm = {Aahad, Jtij, A/^} are included as fit 
parameters and allowed to vary with their current experimental resolutions (t(/sm). We do 
not include as as a fit parameter, which would have only a minor impact on the analysis. 

Formulating the fit in this fashion has the advantage that the probability, P(x^, iVdof), 
properly accounts for the number of degrees of freedom, iVdof , in the fit and thus represents 
a quantitative and meaningful measure for the "goodness-of-fit." In previous studies [1], 
P(x^,7V<iof) has been verified to have a fiat distribution, thus yielding a reliable estimate 
of the confidence level for any particular point in parameter space. All confidence levels 
for selected model parameters are performed by scanning over the desired parameters while 
minimizing the function with respect to all other model parameters. The function values 
where x'^i^) is found to be equal to Xmin + determine the confidence level contour. 
For two-dimensional parameter scans we use Ax^ ~ 2.23(5.99) to determine the 68%(95%) 
confidence level contours. Only experimental constraints arc imposed when deriving confi- 
dence level contours, without any arbitrary or direct constraints placed on model parameters 
themselves. This leads to robust and statistically meaningful estimates of the total 68% and 
95% confidence levels, which may be composed of multiple separated contours. 

The experimental constraints used in our analyses are listed in Table 2 in [8]. It should 
be noted that we use of the e+e~ determination of the SM contribution to {g — 2)^ [23] . 
flSUSY ^ ^3Q_2 ± 8 8) X 10"^°. Indeed, the {g - 2)^ hint has even strengthened with the 
convergence of the previously discrepant SM calculations using low-energy e~^e~ and r decay 
data mmj. 

The numerical evaluation of the frequentist likelihood function using the constraints has 
been performed with the MasterCode [IHTO]. which includes the following theoretical codes. 
For the RGE running of the soft SUSY-breaking parameters, it uses SoftSUSY [25], which 
is combined consistently with the codes used for the various low-energy observables. At the 
electroweak scale we have included various codes: FeynHiggs j26| is used for the evaluation 
of the Higgs masses and a^^^^ (see also PTH^O] ). For fiavor- related observables we use 
SuFla [31] as well as Super Iso ^2]: a-^d for the electroweak precision data we have included 
a code based on [33]. Finally, for dark-matter-related observables, MicrOMEGAs [34] and 
SSARD [35] have been used. We made extensive use of the SUSY Les Houches Accord [36] 
in the combination of the various codes within the MasterCode. 
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The MasterCode framework is such that new observables can easily be incorporated via 
new 'afterburners', as we discuss below for the LHCi/fb constraints. We use a Markov Chain 
Monte Carlo (MCMC) approach to sample the parameter spaces of supersymmetric mod- 
els. The sampling is based on the Metropolis-Hastings algorithm, with a multi-dimensional 
Gaussian distribution as proposal density. The width of this distribution is adjusted during 
the sampling, so as to keep the MCMC acceptance rate between 20% and 40% in order to 
ensure efficient sampling. It should be noted that we do not make use of the sampling density 
to infer the underlying probability distribution. The results of Refs. [8l|9] arc based on a ba- 
sic resampling of the CMSSM with 7 10^ points and a resampling of the NUHMl with 7 10^ 
additional points, both extending up to ?no,mi/2 = 4 TeV. We check that the afterburners 
we apply do not shift the likelihood distributions outside the well-sampled regions. 

3 SUSY searches at the LHC 

The updates concerning SUSY searches in Ref. [8] are based on the public results of searches 
for supersymmetric signals using ~ 1/fb of LHC data analyzed by the ATLAS and CMS 
Collaborations and ^ 0.3/fb of data analyzed by the LHCb Collaboration. For our purposes, 
some of the most important constraints are provided by the ATLAS [16] and CMS [18] 
searches for jets -I- events without leptons, as well as searches for the heavier MSSM 
Higgs bosons, H/A ^7\\19\. Also important are the new upper limits on BR(i?s fi~^^~) 
from the CMS [20], LHCb [21] and CDF Collaborations [37], which arc incorporated in 
Ref. [H]. 

The CMS and ATLAS Collaborations have both announced new exclusions in the (toq: 'm-i/i) 
plane of the CMSSM based on searches for events with jets +1^t unaccompanied by charged 
leptons, assuming tan/3 = 10, = and /i > 0. The updated CMS ar analysis is based 
on 1.1/fb of data [18], and the updated ATLAS 0-lepton analysis is based on 1.04/fb of 
data [16j . It is known that 0-lepton analyses are in general relatively insensitive to the 
tan/9 and Aq parameters of the CMSSM, as has been confirmed specifically for the CMS 
Ot analysis, and they are also insensitive to the amount of Higgs non-universality in the 
NUHMl. Therefore, we treat these analyses as constraints in the {mo,mi/2) planes of the 
CMSSM and NUHMl that are independent of the other model parameters. 

The CMS and ATLAS 0-lepton searches are most powerful in complementary regions 
of the (too7"Zi/2) plane. Along each ray in this plane, we compare the expected CMS and 
ATLAS sensitivities, select the search that has the stronger expected 95% CL limit, and 
apply the constraint imposed by that search 0. This subsequent evaluation/application of 
additional/new contributions is called the 'afterburner'. We assign Ax^ = 5.99, corre- 
sponding to 1.96 effective standard deviations, along the CMS and ATLAS 95% 0-lepton 
exclusion contours in the (jtiq 7^1/2) plane. In the absence of more complete experimental 
information, we approximate the impact of these constraints by assuming that event num- 
bers scale along rays in this plane cx A^^^ where = ^niQ + rn^^^, as described in [7]. 

We then use these numbers to calculate the effective numbers of standard deviations and 
corresponding values of Ax^ at each point in the plane. This procedure has been validated 
by comparing the likelihood it yields with results obtained independently using the generic 



^ It would also facilitate the modelling of LHC constraints on supersymmetry if the results from different 
Collaborations were combined officially, as was done at LEP, is already done for BR(_Bs — > fi~) searches, 
and is planned for Higgs searches. 
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detector simulation code DELPHES, which has also been shown to reproduce quite accurately 
the likelihood function evaluated by the CMS Collaboration using their data [55] . 

4 Searches for a SM-like Higgs boson at the LHC 

Within the supersymmetric frameworks discussed here, a confirmation of the excess reported 
by ATLAS and CMS [22] and consequently the discovery of a SM-like Higgs boson is expected 
to be possible during this year, with a mass in the range between 114 and 130 GeV [22]. 
We assume that this measurement will yield a nominal value of Mh within this range, with 
an experimental error that we estimate as ±1 GeV. In Ref. [5] the possibility that the LHC 
experiments confirm the excess reported around 125 GeV and indeed discover a SM-like 
Higgs boson was analyzed. Assuming 

Mh = 125 ± l(exp.) ± 1.5(thco.) GeV , (2) 

this new constraint was incorporated using the 'afterburner' approach discussed above (see 
also Ref. [8]). 

5 Results from updated SUSY fits 

We now review the effects on the global likelihood functions in various CMSSM and NUHMl 
parameter planes from the inclusion of the latest SUSY and Higgs search data. Below also 
the implications for various physics observables are studied. The (tooj ™i/2) planes shown in 
Fig. [1] are for the CMSSM (left) and NUHMl (right). The regions preferred at the 68% CL 
are outlined in red, and those favoured at the 95% CL are outlined in blue. In the upper 
row the sohd (dotted) lines include (omit) the LHCi/fb data. The open green star denotes 
the pre-LHC best-fit point [5], whereas the solid green star indicates the new best-fit point 
incorporating the LHC SUSY search results. In the lower row the solid (dotted) lines include 
(omit) the assumed LHC Higgs constraint. The open green star denotes the pre-Higgs best- 
fit point [5], whereas the solid green star indicates the new best- fit point incorporating a 
Higgs-boson mass measurement at 125 GeV. 

One can sec in Fig. [Ijthat both new experimental results have a similar, adding up effect. 
The preferred regions arc shifted by both new constraints to substantially larger mi/ 2 and 
somewhat larger toq values. A similar effect can be observed in Fig. [21 where we show the 
(TOi/2,tan/3) planes in the CMSSM (left) and the NUHMl (right), including (omitting) 
the LHCi/fb constraints in the upper row, including (omitting) the the hypothetical LHC 
measurement Mh = 125 ± 1 GeV in the lower row. The preferred values of tan/3 arc shifted 
to substantially higher values. 

The results for the {MA.taaP) planes in the CMSSM and the NUHMl are shown in 
Fig. [3] Again both new experimental results show a similar effect. We observe a strong 
increase in the best-fit value of Ma in both models, especially in the CMSSM, where after the 
inclusion of the latest SUSY searches and the hypothetical Higgs signal now Ma ^ 1600 GeV 
is preferred. However, the likelihood function varies relatively slowly in both models, as 
compared to the pre-LHC fits. 
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Figure 1: The (mo,mi/2) planes in the CMSSM (left) and the NUHMl (right). The Ax^ = 
2.30 and 5.99 contours, commonly interpreted as the boundaries of the 68 and 95% CL 
regions, are indicated in red and blue, respectively. In the upper row the best-fit point after 
incorporation of the LHCi /fb constraints is indicated by a filled green star, and the pre-LHC 
fit by an open star. The solid lines include the LHCi/{\^ data and the dotted lines showing 
the pre-LHC fits. In the lower row the solid lines include the hypothetical LHC measurement 
Mfi = 125±1 GeV and allowing for a theoretical error ±1.5 GeV, and the dotted lines repeat 
the contours including the LHCi/f\^, but without this Mh constraint. Here the open green 
stars denote the pre-Higgs best-fit points whereas the solid green stars indicate the new 
best-fit points. 



6 Implications for the ILC(IOOO) 

In view of the interest in building an e^e~ collider as the next major project at the en- 
ergy frontier, we briefly review the post-LHCi/fb predictions for expectations for sparticle 
production in e+e" annihilation within the CMSSM and NUHMl. In this respect it has to 
be kept in mind that the LHC searches are mainly sensitive to the production of coloured 
particles, whereas lepton colliders will have a high sensitivity in particular for the produc- 
tion of colour-neutral states, such as sleptons, eharginos and neutralinos, as well as yielding 
high-precision measurements that will provide indirect sensitivity to quantum effects of new 
states. 
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Figure 2: The (mi/2,tan/3) planes in the CMSSM (left) and the NUHMl (right), including 
(omitting) the LifCi/tb constraints in the upper row, including (omitting) the the hypothetical 
LHC measurement Mt = 125 ± 1 GeV in the lower row. The notations and significations 
of the contours are the same as in Fig. [IJ 



Fig. 13] compares the likelihood functions for various thresholds in the CMSSM (upper 
panel) and the NUHMl (lower panel), based on the global fits made using the LHCi/fb f^nd 
XENONIOO constraints. The lowest thresholds are those for e+e^ — !■ XiXii Tifi, eaeR and 
P'B.t'-R (the latter is not shown, it is similar to that for eRen). We see that, within the CMSSM 
and NUHMl, it now seems that these thresholds may well lie above 500 GeV, though in 
the CMSSM significant fractions of their likelihood functions still lie below 500 GeV. The 
thresholds for x5x2 s-'^d enCL + glGr are expected to be somewhat higher, possibly a bit 
below 1 TeV. The preferred value for the threshold for lies at about 1700 GeV in both 

the CMSSM and NUHMl scenarios, that for the HA threshold lies above 1 TeV, and that 
for first- and second-generation squark-antisquark pair production lies beyond 2.5 TeV in 
both models. 

What are the impHcations for the ILC(IOOO)? 

Two aspects are important for the correct interpretation of the results reviewed above with 
respect to the ILC(IOOO). First, anything inferred from the coloured sector concerning 
the uncoloured sector depends on the underlying model assumptions, and in particular on 
assumptions about the possible universality of soft supersymmetry breaking at the GUT 
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Figure 3: The {Ma, ta.n 13) planes in the CMSSM (left) and the NUHMl (right), including 
(omitting) the LHCi/f\^ constraints in the upper row, including (omitting) the the hypothetical 
LHC measurement Mh = 125 ± 1 GeV in the lower row. The notations and significations 
of the contours are the same as in Fig. [IJ 



scale. Non- universal models, e.g., low-energy supersymmetric models, or models with dif- 
ferent GUT assumptions, could present very different possibilities. 

Second, while within the CMSSM and the NUHMl the preferred SUSY and heavy Higgs 
mass ranges move to substantially higher values, this upward shift goes along with a sub- 
stantial increase in x^/d.o.f., as can be seen in Tab. [1] Here we review the results for the 
best-fit points in the CMSSM and in the NUHMl "pre-LHC" (i.e. no SUSY searches at the 
LHC, nor any assumption about about any possible Higgs mass measurement), "prc-Higgs, 
post-LHC" (i.e. including the LHCi/fb SUSY searches, but no assumption about any pos- 
sible Higgs measurement) and ~ 125 GcV" (in addition the assumption of Eq. ([2])). 
The droop in the "Fit probability" is clearly visible. The reason is that the pre-LHC data 
favors relatively low SUSY mass scales (to a large extent driven by (g — 2)^, supported by 
M\Y and other observables), while the non-observation of colored SUSY particles favores 
a heavier spectrum, leading to an increasing tension within the CMSSM and the NUHMl. 
Concerning the production thresholds shown in Fig. |4l further increases in the excluded 
regions would yield even higher thresholds, but would also make the CMSSM or NUHMl 
seem even less likely. The time might come to take another look at other GUT based or 
non-minimal supersymmetric models, in which the colored sector (mainly searched for at 
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Figure 4: The x likelihood functions for various pair-production thresholds in e~^e~ , as 
estimated in the CMSSM (upper panel) and the NUHMl (lower panel) after incorporating 
the XENONIOO \39^ and LHCi/fa constraints. The likelihood function for the fiRfiR threshold 
(not shown) is very similar to that for crcr. 



Model 


Minimum 


Fit Prob- 


"^l/2 


mo 


^0 


tan /? 




xVd.oi. 


ability 


(GeV) 


(GeV) 


(GeV) 




CMSSM 














pre-LHC 


21.5/20 


37% 




qo+220 




151^^ 


pre-Higgs, post-LHC 


28.8/22 


15% 


780 


450 


-1110 


41 


Mh ~ 125 GcV 


30.6/23 


13% 


1800 


1080 


860 


48 


NUHMl 














pre-LHC 


20.8/18 


29% 


040+280 




ron+750 


13lf 


pre-Higgs, post-LHC 


26.9/21 


17% 


730 


150 


-910 


41 


Mh ^ 125 GeV 


29.7/22 


13% 


830 


290 


660 


33 



Table 1: Comparison of the best- fit points found in the CMSSM and NUHMl pre- 
Li/Ci/fb JS^, pre-Higgs ^ and including the latest results from SUSY and Higgs searches. 
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the LHC) and the uncolored sector (which is relevant, e.g., for {g — 2)^) are less strongly 
connected as in the CMSSM or NUHMl. 
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